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The gauge invariant, two-component optical conductivity model, with a correlation gap struc- 
ture related to the umklapp scattering processes, is applied to the quasi-one-dimensional electronic 
systems and compared to the recent measurements on the Bechgaard salt (TMTSF)2PF6. The 
optical response of both the insulating and metallic state is found for the half-filled conduction 
band, depending on the ratio between the correlation energy scale 2A2 and the transfer integral in 
the direction perpendicular to the conducting chains, t'^i . The estimated value 2/^%/tf^i agrees rea- 
sonably well with the previous experimental and theoretical conclusions. Parallel to the chains the 
thermally activated conduction electrons in the insulating state are found to exhibit an universal be- 
haviour, accounting for the observed single-particle optical conductivity of the ordered ground state 
of charge-density-wave systems. The band parameters and the related damping energies suitable to 
the normal metallic state of (TMTSF)2PF6 are estimated from the measured spectra. Not only the 
spectral weights but also the damping energies clearly indicate an opening of the correlation gap in 
the charge excitation spectrum. 
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I. INTRODUCTION 

The electrodynamics of quasi-one-dimensional (QID) 
and two-dimensional (2D) electronic systems has been 
the subject of intensive experimental investigation P, |3, 
S H nil 0, II i, 113, in part because various correla- 
tion and dimerization gaps in the charge excitation spec- 
trum are observed recently in the measurements con- 
ducted on several high quahty single crystals. Strong 
evidence for an unusual metallic state in the Bechgaard 
salts (TMTSF)2X, where X = PFg, AsFg or CIO4, has 
been found not only at temperatures close to the critical 
spin-density-wave (SDW) temperature, Tsdw, but also 
at temperatures well above Tsdw Q • 

The anomalies observed in (TMTSF)2PF6 in the latter 
temperature range are the focus of the present analysis 
(the data measured at T « 20 K and T « 100 K 
will be considered here as the typical representatives of 
the T w TsDW and T 3> Tsdw spectra). The most pro- 
nounced anomalous features characterizing the spectra 
measured for the electromagnetic fields polarized in the 
highly conducting direction are as follows, (i) The spec- 
tral weight of the zero-frequency peak is almost indepen- 
dent of temperature, with the related effective number 
of conducting electrons nf^^^^^ ~ n/100 indicating the 
deconfined (metallic) state of the electronic system {n = 
1/Vo is the hole concentration, with Vq being the prim- 
itive cell volume), (ii) There is a pronounced deviation 
of the zero-frequency conductivity from a simple Drude 
behaviour which points at substantial frequency corrcc- 
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tions in the "intraband" damping energy ^iFintra- (iii) 
The maximum in the mid-infrared (MIR) optical con- 
ductivity shifts from ijjpoak/(27rc) = 250 cm"-'^ at T = 20 
K to cjpeak/(27rc) = 800 cm-i at T = 300 K, which is 
approximately the temperature dependence expected for 
the "interband "damping energy ^jFintcr- (iv) This peak 
broadens with temperature exhibiting both a strong sub- 
gap conductivity (with fiFinter ~ 30 meV at room tem- 
peratures, estimated in Ref. ^3|) and a non- universal 
(temperature-dependent) power-law behaviour of the op- 
tical conductivity in the frequency range lo 3> ij-'pcak- 

In order to discuss these anomalies and the existing 
theories describing the electrodynamics of the stro ngly 
correlated strictly ID or QID systems [H El El El iH 
IHElElEHHl, we will consider here a simple, single- 
particle model of the conduction electrons in the presence 
of both two gaps and external electromagnetic fields. It 
is already shown that such a model gives rise the exact 
solution of the site-energy dimerization problem (for ex- 
ample, the anion ordering in (TMTSF)2C104) j2ll, as 
well as a good approximation for the weak bond-energy 
dimerization (and the related single-particle properties of 
the crystals with the charge-density-wave (CDW) ground 
state) E3- Here we will show that this model also pro- 
vides a reasonable description of the correlation effects 
related to the umklapp scattering processes. Unhke the 
exact treatment of the umklapp in the strictly ID ap- 
proaches based on the bosonization procedure [isl , 
suitable to temperatures close to TsdWi the present 
model describes the electrodynamics of the interacting 
electronic system at T 3> Tsdw in terms of the effective 
fermions, and, consequently, is the subject of easy com- 
parison with similar fermionic analyses concerned with 
the room temperature electrodynamic features or 



with other (anomalous) properties of the normal metallic 
state of these salts [1121123.12^. 

The article is organized as follows. In Section 2 and 
Appendix A we present the two-component optical con- 
ductivity for the QID electronic system with the dimer- 
ization/correlation gap in the single-particle excitation 
spectrum, which is a simple extension of the memory- 
function formalism [isl l25l |2^ |27|| . The confinement of 
conduction electrons to the highly conducting chains is 
discussed in Section 3, by considering the spectral weights 
of the zero-frequency and MIR optical conductivity as 
a function of the crossover parameter, i.e. of the ratio 
between the gap magnitude and the transfer integral in 
the direction perpendicular to the chains. The interplay 
among the small parameters of the two-component opti- 
cal model is illustrated then for the insulating and metal- 
lic state of the half-filled band. The results are compared 
to the experimental spectra of (TMTSF)2PF6. The con- 
cluding remarks are given in the last section. 



II. THEORETICAL MODEL 

In this section, the memory-function theory will be ap- 
plied to determine the response of a simple QID metallic 
system to the external electromagnetic fields. The elec- 
trons are assumed to be affected by two single-particle 
perturbations Am(k) with the commensurate wave vec- 
tors Q„i — (27r/(ma), 0), m — 1 and 2 (where Q2 « 2kF). 
a. — ax and b' — b'y are the primitive vectors of the rect- 
angular Bravais lattice, and the z-axis is along the highly 
conducting direction. For the Bechgaard salts, Ai(k) 
is the potential related to the bond-energy dimerization 
along the cc-axis, and A2(k) is an effective perturbation 
which will model the opening of the pseudogap in the 
charge excitation spectrum caused by the umklapp scat- 
tering processes. However, it must be noted that the 
present consideration will be restricted to temperatures 
well above Tgow (^sdw « 12 K in (TMTSF)2PF6 ||), 
and, consequently, the nesting features of the Fermi sur- 
face responsible for the SDW instability of the electronic 
system will be disregarded. 

In order to make the generalization of the memory- 
function approach more straightforward, the analysis is 
divided into three steps. For the crystallographic data 
characterizing the Bechgaard salts, the direct effects of 
Ai(k) on the current-current correlation function are 
found to be negligible 'Tl'|. This issue is briefly recon- 
sidered in Section 2.1. In Section 2.2 the corresponding 
indirect effects are parametrized in terms of the umk- 
lapp interaction gi/2 |28j . and treated by means of the 
mean-field approximation. The current-current correla- 
tion functions are determined then for the case of low 
impurity concentration. Finally, after ensuring that the 
resulting optical conductivity obeys the causality prop- 
erties, the residual scattering processes are taken into 
account in Section 2.3. 




FIG. 1: The hole dispersions ei{kx,n / {2b')) (solid lines) as a 
function of the perturbations Ai(k) (a) and A2(k) (b). The 
Ai(k) 7^ dispersion is contrasted to the weak-dimerization 
(WD) and strong-dimerization (SD) dispersions defined in 
the text. The dot-dashed line is the Fermi level in the 
quarter-filled Ai(k) = band, and £f = eo(7i"/(2a), 7r/(2b')) — 
eo(0, 7r/(26')) is the related Fermi energy. The dispersions of 
the A2(k) = and A2(k) 7^ SD models are shown in the 
reduced zone representation; c, c and C, C_ are the corre- 
sponding band indices. 



A. Bond-energy dimerization 

The bond-energy dimerization along the highly con- 
ducting chain axis occurs in all Bechgaard salts, and is 
usually described in terms of the dimensionless parameter 
a = (ai — a2)/a, with ai and 02 being nearest-neighbor 
distances and a = ai + 02. In (TMTSF)2PF6 the char- 
acteristics of the triclinic Bravais lattice are found to be 
ai-a2 = 0.03 A, a = 7.297 Kb' = bsin-f = 7.291 A and 
a - (b X c) = 714.3 A^ 

In these compounds the conduction electrons are usu- 
ally described using the quarter-filled hole band eo(k) = 
— 2iaCosik • a — 2tb'C0sk • b', with the perturbation 
Ai(k) = i(7asinik • a. Here and th' are the bond- 
energies for a — 0, and g is the related electron-phonon 
coupling. After diagonalization the lower band becomes 
half-filled, with the dispersion e(k) = (l/2)[eo(k)-|-eo(k± 
Qi)] - V(l/4)[eo(k) - eo(k ± Qi)]^ + |Ai(k)|2, as iUus- 
trated in Fig. 1(a). 

The spectral weight of the interband optical excita- 
tions relative to the total spectral weight is given ap- 
proximately by [A]'/|e(kF)|]^, with Aj* = ga representing 
the magnitude of the dimerization potential and Egi « 
|2e(kF)| is the related threshold energy (see the differ- 
ence between the total and intraband spectral weights 
for the case A(k) Ai(k) in Eqs. lEHll, lEU, (E3) and 
(IX6|)). At least in (TMTSF)2PF6 this ratio is negligible 
[Il|. Therefore, the direct effects of Ai(k) on the opti- 
cal conductivity will not affect the forthcoming analysis 
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and we can write the bare hole dispersion in the form 
e(k) « eo(k) (hereafter, this limit is referred to as the 
weak-dimerization (WD) limit, see Fig. 1(a)). Conse- 
quently, the quantitative analysis of the spectral weights 
in Sec. 2.3 will employ the WD dispersion. On the 
other hand, for the qualitative consideration of the low- 
frequency spectra, the strong-dimerization (SD) disper- 
sion is more adequate due to simple symmetry properties 
of both £(k) and the related current and Raman vertex 
functions. 



B. Umklapp scattering processes in the mean- field 
approximation 



Hq is the bare Hamiltonian. H[ represents the single- 
particle scattering processes (due to impurities, phonons, 
CDW or SDW fluctuations) and H'^^^ is the coupling 
Hamiltonian which couples the electromagnetic fields to 
the current and Raman density fluctuations. In the 
memory-function approach at temperatures well above 
TsDW the residual two-particle scattering processes, rep- 
resented here by H!^^ play usually the marginal role and 
therefore will be neglected in the present section. 

The bare Hamiltonian comprises two terms; the first 
one (7^0 ) includes the contributions which are diagonal 
in the Zk representation, and the second one (Hq) the 
off-diagonal contributions coming from the single-particle 
potential A(k) = e*'^«|A(k)|: 



From the work of Emery et al. [23 it is known that 
the electrodynamic features of the Bechgaard salts are 
affected by the bond-energy dimerization dominantly 
through the umklapp interaction gi/2 oc UAi/ep. Their 
conclusions have been supported by the systematic inves- 
tigation of the optical conductivity in the MIR part of 
the spectra j^. It has been established that, as a shifts 
to higher values, the energy i?poak — ^peak becomes 
larger going from Speak « 25 meV in (TMTSF)2PF6 
{a « 0.004 i2a]) to Speak « 200 meV in (TMTTFjzPFg 
(a « 0.014)). Additional evidence for the crucial role 
of these scattering processes was given by Favand and 
Mila [l3|, considering the MIR optical features of the 
dimerized ID Hubbard model. They have found that the 
dimerization potential, together with the strong Hubbard 
interaction U, results in a renormalized threshold energy 
Eg2 ~ 2A]', in addition to the bare threshold energy Egi 

{Eg2 < -Bgi). 

In order to study the influence of these processes on the 
electrodynamics of the conduction electrons described by 
the (WD or SD) hole dispersion £(k) at T > Tsdw, 
we adopt now the mean-field approximation, with the 
effective perturbation A2(k) (x 51/2- 



1. Model Hamiltonian 

As illustrated in Fig. 1(b), the electrons are shown in 
the zone representation \kx\ < 7r/(2a), \ky\ < ir/b' (for 
simphcity, the 2D Brillouin zone of (TMTSF)2PF6 is ap- 
proximated by the rectangular zone). The band indices 
of A^ = and A^ 7^ wiU be denoted by I & {c,c} 
and L € {C,C}, respectively. The ei(k) and i?L(k) are 
the related hole dispersions (obviously, ec(k) = e(k) and 
ec(k) = £(k ± Q2)). For the sake of generality, the index 



TO = 2 (in A2(k), Q2 and Eg2 = Ecikp) - Sc(kF); Eg2 
is the MIR threshold energy) will be omitted and an arbi- 
trary filling of the conduction band will be assumed. The 
final results (Section 3) will be given for Q Q2 w 2kF, 
with the index to = 2 explicitly written 
The total Hamiltonian reads as 



H = Ho + H[ + H!^ + H"" 



(1) 



Ho = Y.^o'mljL, 



(2) 



HS%k) = ee(k),^ iff (k) = £,(k) and i7o^^(k) = A(k). 
The transformations 



^!7k(/,L)LL, 



lead to 



Ho - 5]SL(k)4^Lk., 



(3) 



(4) 



Lk(T 



with the dispersions 



i?c,c(k) = i[ee(k)+e.(k)]± Jie2^(k)-f |A(k)p. 



(5) 



The transformation-matrix elements are given by 



/C/k(c,C) Ukic,C) 
I [/k(c,C) (7k(c,C) 



^(k) 
cos — - — 

-e'^osin^ 
2 



^(k) 



cos ■ 



^(k) 



(6) 



where 



tan V3(k) 



2|A(k)| 

ecc(k) 



(7) 



and ecc(k) = ec(k) - £c(k). 

In the Lk representation, one can write 

H[ ^ J2 ^"'(k,k')LLikv, 



(8) 



LL'kkV 



H° 



E 

,L'k, 

x^^(q)7r(k;2) 



LL'ka 

l2 f„\^.LL 



c Zmc 



(9) 
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The electromagnetic fields are described in terms of the 
transverse vector potential Aa(r) polarized in the a di- 
rection {Aa{q) and A^(q) are the Fourier transforms 
of ^^(r) and ^^(r), respectively), while J^^ (k) and 
7^^(k;2) are the related current and bare Raman ver- 
tex functions (see Appendix A). 

The memory-function approach will be applied now to 
the Hamiltonian (Q. 



2. Current-current correlation functions 

The retarded current-current correlation functions are 
defined by H 



1 



with ,/^^ (q) representing the current density operator 



tL'Lf 



Ji'^'(k,k+) = Ji:^'(k)iLW- 



LL' I 



(11) 



Here k+ is the abbreviation for \<i + (\, z — hw + vq and 
{{A\ B))z labels the retarded correlation function of the 
density operators A and B: 

/oc 
dt e-'^"\{A;B))u 

{{A-B))t = -m){[m.B]) (12) 

[Ait) is the Heisenberg representation of A). 

In order to determine 11^^ (q, z), which is now a rather 
complicated function of various scattering terms (ffp, H[ 
or iJj), one usually starts with the equation of motion 
for the correlation function 



(13) 



I?f^'(k,k',q,<) = ((4"(k,k+);ji'^(kV,k')))* 
which gives 

ZLL,(k)I?f^'(k,k',q,z) (14) 



= <5k,k'?i|J„"'(k)|2[/4k) - /L,(k+)] - 



1 



^LL'(k') 



x[P2" (k,k',q, Zii,(k'))-I?2" (k,k',q,0)]. 

The first term on the right-hand side of this expression 
is relevant only for the interband processes, and con- 
sequently was absent in the previous memory-function 
analyses [H IM El 113 ■ In this term the effects of H'^ 
on the photon absorption/emission are described using 
the exact diagonalization of the single-particle problem, 
Eqs. ©-Q, with the momentum conservation in the in- 
terband current vertices fulfiled due to Q in Hq. On the 



other hand, the second term defines the intraband re- 
laxation function, treating the relaxation processes (pro- 
portional to {H')^, H' ~ -Hj) in the way similar to the 
Boltzmann equations [sj, |33| • 

Here T>2^ (k, k.' ,q, ZTj/Ck')) is the force-force correla- 
tion function [23, l27l associated with the current- 
current correlation fimction Vf^ (k, k', q, z): 

I?2"'(k,k',q,ZLL,(k')) (15) 

= (([ji'^'(k,k+),i?'];[ji''^(kV,k'),i?'])).. 

Furthermore, ZLL'(k.) = z + i?L(k) — £'L'(k+) and 
/L(k) = [I + e/'I^^^k)-,.]]-! -g ^jjg Fermi-Dirac distri- 
bution function. 



3. High-frequency limit 

It is rather a good approximation to retain in the inter- 
band relaxation function only self-energy contributions 
[2l|| . show this term in the form — iyyintor^^f ^ (k, k', q, z), 
neglect in Eq. 1)141) the terms which are proportional 
to q and, finally, focus attention on the simplest form 
of the scattering Hamiltonian where the interactions 
l/'-^'-^(k, k') are time- independent. 

The leading contributions in H' to ^ (k, k', q, z) are 
given now by 

zl?f^(k,k',q,z) = -i[l?C^^(k,k',q,z) 

z 

-I?2^^(k,k',q,0)], 
Pf^(k,k',q,z) = <5k,k'kf^(k)|'l?cc(k,k,z), 

(16) 

with the high-frequency intra- and interband current- 
current correlation functions {LL' = CC and LL' = 
CC_, CC, respectively) of the form 

zli-T'-'iz) = ^Y. ^^^^^^^^[Jf^(k) 

kkV ^ 

- J„^^(k')] 'i[I?cc(k, k', z) - VcciK k', 0)], (17) 
niir'°°(^) = ^E['^"''-(k)]'^fc(k,k,z) 

ko- 

+I?cc(k,k,z)] (18) 

(the label oo stands for the high-frequency limit 
^/'7intra ^ 1; ?7intra IS the intraband damping energy to 
be defined below). All physically relevant contributions 
are shown here in terms of 



I?ii,(k,k',z) = ((iLikv;4.ik.)). 

, /L(k)-/i,(k') 



z + EL{\^)~EL'{k'y 



(19) 



For example, in the single- impurity-scattering approxi- 
mation, one obtains (|T/'^'^(k, k')H « A^,/iV|F°(k, k')P; 
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Ni is the presumably small number of impurities and 
y''(k, k') is the scattering potential due to one impurity 
[3^ . The generalization of Eq. I|17|l to the case where 
the interactions F'-^'-^(k, k') are time-dependent or where 
the two-particle scattering term H2 is taken into account 
explicitly is straightforward and will not be considered 
here (see, for example, Ref. 



4- Low-frequency limit 

In the memory-function approach to the normal metal- 
lic state, the low-frequency current-current correlation 
functions come on collecting the most singular scattering 
events in powers of {H'Y /lo. The standard procedure 
|2,'t| is based on the intraband memory (or relaxation) 
function which is defined by 



Ma (2) 



intra, 00 



(^) 



/i-i cff 
"intra. a 



(20) 



,eff 



is the effective number of conduction electrons 



^eff 



^2ylj-^^ ( (k) 

kcr 



(21) 



kcr 



with 7£^(k) being the static Raman vertex, Eq. HA6|I . 

In the metallic systems with the anisotropic electron 
dispersion, the numerator of Eq. (|20ll has a complicated 
structure. Nevertheless, for ?kj — > 0, a formal decou- 
pling of two integrals in the expression (|17() is possible if 
the average-inverse-relaxation-time approximation is ap- 
plied. In this case, the inverse relaxation time 



E(l^^^(k,k')P)[l~jf'^(k')/J.^^(k)] 



'nntra(k) 



x27r(5[£;c(k') -Ai] 



(22) 



is averaged over the Fermi surface, resulting in the intra- 
band damping energy 



»7intra = Im{ Af^ (0 + Ita) } = 



(k) 



). (23) 



Finally, the resulting intra- and interband current- 
current correlation functions read as 



n" 



„2 cff 
"intra. a 



l^intra 



m ftw + i77intra ' 

nL"r'°°(fi^ + iw)- (24) 



A more rigorous treatment of the intraband relaxation 
processes should include the explicit numerical calcula- 
tion of two coupled integrals in the expression H17|l : how- 
ever, this is beyond the scope of the present work. 



as the other relaxation processes not considered in the 
?7i ^ model) 0, 113 • Depending on the ratio 



C. Optical conductivity 

Due to causality, the structure of the multi-component 
optical conductivity comes on combining the expressions 
for the current-current correlation functions obtained for 
77i ^ with the corresponding diamagnetic contribu- 
tions, and replacing r\i with the frequency-independent 
damping energies ^iFi (representing the single-impurity 
relaxations for the usual impurity concentrations as well 
at ion pr o 

among the band parameters, this procedure will result 
in different optical-conductivity expressions. The two- 
component optical conductivity (where 2A°/£'g « 1) and 
the generalized Drude formula (2A°/ii^g <C 1) are the 
most interesting cases. 

As a first approximation it seems natural to interpret 
the optical conductivity of (TMTSF)2PF6 in terms of 
the two-component model with the correlation energy 
scale 2A° comparable to the threshold energy and 
with both ?iFintra and ?iFi,itcr independent of frequency. 
The total optical conductivity CTjj°''*'(a;) is the sum of the 
intra- and interband contributions, which are given by 
the gauge-invariant expressions 0, 0, l3^ 



^intra/ 
inter/ \ 



■ p2 efT 

1 ^ "intra. a 



LU m huj + i/lFintra 
i 1^ (M2|J^C(k)|2 



(25) 



ujV 



/c(k)-/c(k) 

huJ - EccO^) + i^rintcr 

/c(k) - /c(k) 
hu + -Ecc(k) + ifi-Fint, 



(26) 



Since the residual scattering processes (corresponding 
primarily to the frequency-dependent contributions in 
H[, or to the two-particle term H2) become increasingly 
pronounced at low frequencies, the intraband contribu- 
tion 125|l should in principle be replaced by the general- 
ized Drude formula 



hLu{m{u!)/m) + ilTa{Ma{hLV + irj)} 



.(27) 



The total optical conductivity is now the sum of the 
expressions (|27ll and H2t)|l . Here m{uj) results from the 
Kramers-Kroning relations, with H' = H[ + H'2 in Eq. 
(Pn|l . Similarly, in the limit 2A°/Eg < 1, one obtains the 
generalized Drude expression crjj°*'''(w) « a™^'^^{uj) with 
H' ^ H'q + H[+ H'2 and CTjf*^'^(w) given by Eq. 

It should be noted that the direct integration of the 
real part of the expressions (|25|l and (|26() gives rise to 
the intraband, interband and total spectral weights of 
the form ITllsl 



1 



m 1 2 

0'"j,a -i^i'O; 

maa 2 



-Van. 



(28) 
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8 (hole doping) 

FIG. 2: The development of the effective number of con- 
duction electrons with doping for the A" = WD and SD 
models. The dot-dashed line is the prediction of the free-hole 
approximation (FHA). = 0.25 eV, tb' — 25 meV (WD case) 
and ta = 0.125 eV, tb' = 12.5 meV (SD case). For clarity, the 
a = b' effective numbers are multiplied by a factor of 10. 

where i G {intra, inter, total}. The effective total number 
of electrons is given by 

<oVa = -^E7of(k;2)/c(k), (29) 

kcr 

'^fnter,a = <ftaLa " "?ntra,a ^ nd nfg^^^^ C OmCS frOm 

Eq. jni- Furthermore, f^o = y/ine^ /{mVo) IS an aux- 
iliary frequency scale which represents the bare plasma 
frequency attributed to the free holes with the concentra- 
tion n = I/Vq, rriaa = 2ft^/(taa^) is the mass parameter 
and the (maa/™)Von^^ are the effective numbers shown 
in a convenient dimensionless form. 

Before turning to the detailed numerical calculation, 
let us consider the total spectral weights and remember 
the estimates of the bond-energies in (TMTSF)2PF6 H 
It has been established that the measured plasma 
energy /i^ltotai.a (1.1-1.4 eV Q, HH is close to the 
energy « 1.385 eV (or 17n/(27rc) « 1.1 • W cm'i). 
Using the A° = WD model |3] this observation has led 
to the conclusion that the bond-energy ~ t^, where 
ia ~ 0.285 eV is the characteristic bond-energy for which 
the mass parameter maa is equal to the bare electron/hole 
mass m. Similarly, for d ~ Ijcorresponding to the open 
Fermi surface), one obtains 

'^totaUb' ~ /3"^'^total,a> (30) 

with /? 4 for the A° = WD case (see Fig. 2). For 
tb' = O.lta, this has resulted in the bare plasma frequency 
^^totai.b' ~ 0.2r2toiaiA which agrees reasonably well with 
measured values H,y,Q. 



Due to the scattering term Hq, for Q « 2kF the 
quantities n?^^^^ and W^^'^'°°{z) are strongly reduced 
with respect to their values in the free-hole approx- 
imation. In addition, the analytical decoupling of 
two integrals in 11'^^'^^'°° (z) can be achieved only us- 
ing the average-inverse-relaxation-time approximation 
{~hLjIm{U^'J-'°°{fiw + ir,)} « hTi^t,,e^n'?g^^ Jm). For 
this reason, from here on we will focus the attention 
on the two-component model H25|) - (|26|l only, with the 
damping energies independent of frequency and with 
A(k) = AO. 

Leaving the factors m/maa and ri§/2 in Eq. H28(l off, 
the total spectral weight is calculated now as a function 
of the doping level for the bond-energies ia — 0.25 eV, 
tb' = 25 meV (WD case) and L = 0.125 eV, ib' = 12.5 
meV (SD case) and for the photon polarizations a — a 
and a = b' , and shown in Fig. 2. Although the SD model, 
when fitted to the measured total spectral weights, will 
result in the slightly different values of the bond-energies, 
on the qualitative side (for not too large A°/ (2fa) and 
huj/{2ia)) there is no essential difference between these 
two models in explaining the structure of the A° 
optical conductivity. In this respect, we continue the 
analysis considering the SD model, and assuming that 
/iw, fceT, A°, tb', fiFi <C 2ta. It is important to notice 
that the model (|25|l - H2t)|l for tb' — > is essentially the 
same as the strictly ID spinless-fermion optical model of 
Pedron et al. j^, and, not surprisingly, the present es- 
timates of ta, A° and HTi will be close to their estimated 
values. However, there is a pronounced qualitative dif- 
ference between these two optical models regarding the 
Fermi energy and wave vector positions (Q « 2kF here, 
contrary to Q w 4kF in the spinless-fermion model). 



III. COMPARISON WITH EXPERIMENTS 

The optical anomalies of (TMTSF)2PF6 mentioned in 
Section 1 will be analyzed now in the framework of the 
QID fermionic model (|25|l - (|26[l . The competition be- 
tween the small energy scales of this model, which re- 
sults in the structures in the optical conductivity simi- 
lar to that observed in experiments, is the focus of the 
present discussion. We consider first a few general fea- 
tures which emerge as a result of this competition. Then 
the model parameters suitable to the T = 100 K spectra 
of (TMTSF)2PF6 wiU be discussed in more detail. In the 
rest of the text the full notation is used with the index 
TO = 2, in A2 and Eg2, written explicitly. 

A. Competition between tb' and A2 

The a — a conductivity sum rule of the present two- 
component model, f^^^t^i^ = l^Ltra,a + ^^fntor,a' ^OCS UOt 

depend on the damping energies hVi. The main conse- 
quences of the interplay between tb' and A2 can be thus 
seen from Fig. 3, where the effective numbers nf^^j.^ ^ 
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FIG. 3: Main figure: The a = a intraband spectral weight as 
a function of the chemical potential /j. for different values of 
2ib/ (2fb' = and 5 - 11 meV, 2^ = 0.25 eV, A§ = 10 meV 
and T = 10 K). Inset of figure: The dependence of nf^^j^ 
on S for the = (dashed line), = 10 meV (solid line) 
and A§ = 20 meV (dotted line) SD models, 2^ = 0.25 eV, 
2tb' = 25 meV and T = 10 K. The a ^ b' spectral weight is 
again multiplied by 10. 



are displayed as a function of the crossover parameter 
2A2/tb' at a relatively small temperature T = 10 K. For 
tb' — > but A" finite, the electrons will be confined on in- 
dividual chains, as found in the insulating (TMTTF)2v'C 
salts, where X = PFg or Br j^. The effective number 
'^fntra a ~ n/100 required to give the metallic behaviour 
observed in (TMTSF)2PF6 is obtained for 2A^/tb' ~ 5 
(main figure, T = 10 K). This value is approximately 
twice the critical value Eg2/th' = 1.8-2.3 for the con- 
finement (i.e. for the metal-to-insulator phase transi- 
tion) obtained by the renormalization group approach 
for the problem of two coupled chains (l^l- The value 
of 2A2/ib' decreases with decreasing temperature, re- 
sulting in 2A^/tb' w 4 at T = 2 K. Another reason 
for this discrepancy should be the fact that the ratio 
tg_/t\ji plays important role as well; namely, for ta fixed, 
the increase of transforms the QID electronic system 
into an anisotropic 2D system, where the dimerization 
along the x-axis becomes less important. As illustrated 
in the inset of the figure, for 2A2/tb' small enough, the 
gap opens only on a part of the Fermi surface result- 



ing for n\ 



eff 



in a local maximum at (5 = 1. On the 



contrary, for 2A2/tb' > 2.5 and T = 10 K, the correla- 
tion gap develops on the complete Fermi surface, with 
the bare group velocity Jf'^(k)/e in the effective number 
of conduction electrons replaced by the "renormalized" 
one, J^'-^(k)/e = cosip{k)Jl'^{k)/e, with approximately 
cosi^(k) « cos93(kF). 



the same way as nf^^^^ ^. While the rest of the a = a to- 
tal spectral weight, i.e. the difference between the dashed 
and solid (dotted) line in the inset of Fig. 3, reappears 
through the coherent excitations across the correlation 
gap £'g2, for a — b' there is no counterpart of these co- 
herent finite-frequency excitations. Therefore, to fulfil 
the a = h' conductivity sum rule, one has to go beyond 
the optical model (|25|) - (|2(j|l . and to take also incoherent 
optical processes in the b' direction into account [39| . 
Obviously, this issue is of a general importance for the 
complete quantitative description of the optical conduc- 
tivity of the Bechgaard salts, and should be in relation 
with the disappearance of the plasma edge in the a — b' 
reflectivity spectra of the (TMTTF)2X salts 0. 



B. Relation with the CDW systems and the 
high- Tc cuprates 

Besides the dramatic dependence on the crossover pa- 
rameter 2A2/ib', the optical model H25|l - H26|) exhibits 
also various temperature-dependent features, due to 
the competition between the energies fceT, A2(r) and 
printer (T")- Some of these temperature effects are at first 
sight similar to the anomalies (i)-(iv) mentioned in Sec- 
tion 1. Interestingly, they are more directly related to 
different MIR optical features observed in various other 
QID and 2D systems. To compare briefly the predictions 
of the model with the MIR spectra measured in the CDW 
systems in the ordered ground state , as well as in the 
normal metallic state of the underdoped high- Tc cuprates 
we apply it now to the insulating (5 = 1, tb' 

case. 

The real part of the intraband conductivity has the 
usual form Re{(Ti'^*-(t^)} = a^^/{l + (c^/Fintra)'). 



(jDC = nlVorifg^^ JiAnTintr.) is the DC conductivity, 
which vanishes at T = OK. The temperature depen- 
dence of (jjp^ is governed by the thermal activation of 
the conduction electrons, according to the expression 



,eff 



-E 



7fa'^(k)[/c(k)-/c(k)]. (31) 



This temperature dependence is accompanied by the de- 
crease of the interband spectral weight, as seen from 



Re{ar'(^)} 



--E 

ujV ^ 



{hiu)^\j^'^{k)\ 



AhwEcc (k) 



printer [/c(k) - /c(k)] 



(32) 



The number 



intra, b' 



, is affected by the correlations in 



where the same temperature factor can be recognized 
(/c(k) - /c(k) ~ tanh[i?c(k)/(2fcBr)] for ^l = 0). 

It should be interesting first to remember that the 
optical model H32|) with A2 — > Acdw gives a correct 




FIG. 4: The dependence of the real part of the a = a optical 
conductivity on k^T, for 5 = 1, 2t-^ = 0.25 eV, t^i —> and 
'irintra ~ 5 mcV. Main figure: The temperature dependence 
of the spectra for A2 Acdw < 10 meV and ftFintcr = 5 
meV (underdamped regime), cr^'^ is the DC conductivity of 
the metallic A2 = state, and the solid line is the related 
optical conductivity. Inset of figure: The dependence of the 
total and interband spectra on fiFintor: Acdw = 4 meV, T = 
56 K, ftFintcr = 5 meV (curves A) and ftFintcr = 2 meV (curves 
B). 



FIG. 5: The dependence of the real part of the a = a optical 
conductivity on the ratio ftFintcr/A2, for S = 1, 2fa = 0.25 
eV, tb' ~* 0, ftFintra = 5 meV, T = 10 K, nis^^/m = 1 and 
ftni/(47r eV) = 0.355 • 10^(J7cm)-\ Main figure: The over- 
damped regime; fiFintcr = 30 meV. The data measured in 
(TMTSF)2PF6 at r = 20 K are given for comparison. In- 
set of figure: The intermediate regime; A2 = 10 meV with 
ftFintcr = 30, 20 and 10 meV for the curves A, B and C, re- 
spectively. The energies ftFintcr (main figure) and 2A2 (inset 
of figure) are indicated by arrows. 



prediction for the T « K subgap conductivity char- 
acterizing the ordered CDW systems, due to the gauge- 
invariance factor (?ia;)^/£'^p(k) [l^. It also provides the 
correct treatment of the conductivity sum rules, even for 
the relatively large damping energies, due to the factor 

4awi;cc(k)/[(fi^ + %c(k))i+/ir2 er]' at variance with 
the ordinary CDW model which fails in both 

of these cases. 

Fig. 4 illustrates the development of Re{(T^°*'^'(ci;)} 
with temperature in the CDW ground state for the un- 
derdamped case Ao/(firintcr) = 2. The temperature de- 
pendence of the order parameter is modeled by a simple 
expression Acdw(T') Aqv^I - T/T* (T* is an auxil- 
iary temperature scale given by ksT* = Ao/1.75, with 
T* « TcDw)- Interestingly, almost the same tempera- 
ture behaviour is found for the a — b' optical conductiv- 
ity in (TMTSF)2PF6 at T < T* w 15 K, with the SDW 
(pseudo)gap parameter 2Ao/(/ic) w 70 cm^^ 6]. 

Equally important is the observation in the normal 
metallic state of Bi2Sr2Cu06 that the MIR maximum 
shifts with increasing temperature from a;poak/(27rc) ~ 
100 cm-i at T = 30 K to Wpoak/(27rc) « 500 cm^^ 
at T = 300 K 9j (with a similar trend found in 
the La2-a;Srj;Cu04 compounds |lOj). This unusual be- 
haviour has been explained using the ID orbital Kondo 
model of Emery and Kivelson [16| . which predicts 
^'peak ~ sr. Fig. 5 demonstrates that the same be- 



haviour emerges in the overdamped regime of the model 

ftFintor > leads to Spoak ~ ^Fintor- It 

seems that at least in the compounds mentioned above 
the unusual temperature dependence of iJpeak at tem- 
peratures up to room temperature has the same origin, 
the competition between the correlation energy scale in 
question and the presumably large (interband) damping 
energy. 

In conclusion, in the underdamped regime of the 
present optical model the position of the maximum in 
the MIR spectra is close to the threshold energy Eg2 as 
well as to the correlation energy scale 2A2. On the other 
hand, in the intermediate and overdamped regime one 
expects to be i^peakCAj, fi.rinter), so that the correlation 
energy 2A2 cannot be extracted directly from the mea- 
sured Speak- 



C. (TMTSF)2PF6 

Once the bond-energies ia and together with the 
correlation energy scale 2A2, are estimated from the 
spectral weights, the damping energies ftFintra and /iFjntor 
and the background dielectric function £00, a are the only 
free parameters to be estimated from the real part of the 
optical conductivity and the real part of the dielectric 
function obtained in experiments. 
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FIG. 6: The real part of the a — a optical conductivity 
as a function of the intraband damping energy, for the SD 
model with nfftra.a = n/100; 2^ = 0.25 eV, 2^, = 8 meV, 
A§ = 10 meV, ftFintor = 30 meV, T = 10 K, m^/m = 1 and 
fing/(47r eV) = 0.355 • 10^(ncm)~^ fiFintra = 1 and 0.1 meV 
for the curves A and C, respectively. The prediction of the 
simple Drude model (curve B, A2 = and fiFintra = 0.1 eV) 
and the typical measured spectra 0| are also shown. 



FIG. 7: The dependence of the real part of the a = a di- 
electric function on ftFintcr (main figure) and /iFintra (inset 
of figure), for 2^ = 0.25 eV, 2ty,, = 8 meV A§ = 10 meV, 
T = 10 K, maa/m = 1, hQ,o = 1.4 eV and e^,a = 2.5. Main 
figure: fiFintra = 1 meV, and fiFinter = 15, 20 and 25 meV 
for the curves A, B and C, respectively. The dot-dashed line 
describes the A2 — case. Inset of figure: ftFinter = 20 meV 
and fiFintra = 0.1, 1, 2 and 5 meV (D, B, E and F). The 
experimental data are from Ref. 2l- 



According to the general relation 

ea(t^) ~ £oo.a(c^) + — ^r"'(c^), (33) 
UJ 

eoo,a(w) — 1 describes the contribution of the optical 
processes not included in the Hamiltonian In the 
Bechgaard salts, for hiu < \ eV, it is established that 
Ini{eoo,a(t^)} ~ and Re{eoc.,a(w)} = eoo,a ~ 2.5 1, 2J. 

In order to obtain a satisfactory explanation of the 
estimated damping energies hVi, it is useful first to re- 
member psl liol l4l| that the equation of motion H14|) 
generates automatically the self-energy and vertex con- 
tributions to the memory functions (i.e. to hVi). Further- 
more, it should be noticed that omitting the vertex cor- 
rections makes the damping energies fiFintra and SFintor 
appear as two limiting cases of a frequency-dependent re- 
laxation function hV{huj) which describes the pseudogap 
effects in the electron self-energy (/iFintra ~ ?iF(0) and 
printer ~ hT{hw 3> SAj)). Not surprisingly, the function 
hT{huj) related to the umklapp scattering processes be- 
yond the mean-field approximation adopted here should 
have the temperature and frequency dependence similar 
to that found recently for the half-filled spinless Holstein 
model in which the CDW (pseudo)gap features have been 
studied 

The dependence of Re{cr^°*'''('^)} and Re{4°*''Hw)} on 
the damping energies HVi for the SD case with nf^^j^ ^ = 
n/100 is shown in Figs. 6 and 7, and compared to the 
experimental data. According to the questions raised in 



Section 1, the most important qualitative conclusions are 
as follows. 

The power-law frequency dependence Re{(T^°*'''(a;)} oc 
w^", V w 1.3, for frequencies greater than £^peak/fij ob- 
served at T < 20 K, is well understood as a direct con- 
sequence of the umklapp scattering processes 0, |^ . 
On the other hand, the temperature dependence of both 
u and E-pcaki in particular at temperatures above 100 
K, indicates the regime where ftFjntcr is comparable with 
2A2, and, consequently, where the temperature effects 
are primarily due to fi.Fintor(^r), as shown in the inset of 
Fig. 5. Therefore, the energy ftFintcr ~ 25-30 meV rep- 
resents roughly the high-frequency limit of the relaxation 
function hT(hhj). On the other hand, the zero-frequency 
limit reveals the energy fi-Fintra ~ 1 meV at T = 100 
K, indicating the strong pseudogap effects, with an ad- 
ditional, pronounced decrease on decreasing temperature 
(IFintra ~ 0.1 nicV at T = 20 K). To improve the overall 
agreement between the model predictions and the mea- 
sured spectra, including the anomalous frequency depen- 
dence of the Drude peak as well as the temperature de- 
pendence of its spectral weight, it is necessary to take 
into account the frequency dependence of the damping 
energies in Eqs. (|25|I -H26 () . When considering the two- 
component model with the causality properties treated 
correctly, this problem, however, requires a more accu- 
rate and extensive investigation and will not be discussed 
here. 
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According to Fig. 7, there are two zeros of Iie{ea{i^)}, 
representing the motion of two different plasma modes. 
These frequencies are related to the frequencies fi^^a es- 
timated from the spectral weights. Not surprisingly, 
a simple Drude relation ujp^^ « ^^i,a/^£oo,a is not 
valid here. In this respect, the total plasma frequency 
/{2ttc) ~ 0.55 • 10* cm~^ has to be compared to 



.total 
pl,a 



the plasma frequency il,o/(2TrCy/el 



0.7 • 10* cm-i 

of the Aq = case, and to the bare plasma frequency 

0, 



1.1 • 10* cm ^ corresponding to 



17o/(27rc) 

£oo,a = 1- Similarly, uj™^1''/{2ttc) « 20 cm"\ in con- 
trast to riintia,a/(27rc) « 10'^ cm^*. Furthermore, a pro- 
nounced dependence of t^pfl'^ on ^iFintra for Fintra > t^pi'i^ 
can also be noticed in the inset of figure. 



usual substitution for the electron momentum, and, 
finally, transform the obtained coupling Hamiltonian in 
the Lk representation. 
Inserting 



SHll'ik) 



dH{l (k) e 
dka he 



^a(q) 



(Al) 



ff-t = J2 S< (k) [li+^jL + h.c] , (A2) 



IV. CONCLUSION 

To summarize, we have formulated a gauge invariant 
approach to the optical conductivity of the QID inter- 
acting electronic systems with a MIR structure clearly 
distinguished from both the zero-frequency (Drude-like) 
and the high-frequency (background) contributions. By 
using the simplest limit of this model, with the frequency- 
independent damping energies, the confinement of the 
conduction electrons to the highly conducting chains is 
analyzed. The estimated ratio 2A2/ib' required to give 
the spectral weights measured in (TMTSF)2PF6 agrees 
with the conclusions of the previous renormalization- 
group study. 

The real part of the optical conductivity is calculated 
for the typical insulating and metallic state, giving rise 
a satisfactory agreement with the experimental findings 
at temperatures well above Tsdw- The estimated damp- 
ing energies clearly show the (pseudo)gap features in the 
charge excitation spectrum and are in a qualitative agree- 
ment with the relaxation function calculated recently for 
the systems with the CDW ground state. However, to ob- 
tain a better fit to the measured data, even at T 3> Tsdw, 
one has to go beyond the present model and to treat the 
umklapp scattering processes dynamically. 

Two plasma modes are found in the metallic state, 
with a nontrivial relation between their frequencies and 
the corresponding spectral weights. 
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APPENDIX A: VERTEX FUNCTIONS 

To account for the effects of the dimeriza- 
tion/correlation potential on the coupling between 
the tight-binding electronic system and the electromag- 
netic fields, we choose the /k representation, apply the 



we obtain the coupling Hamiltonian given by Eq. ^ in 
the main text. The intraband bare Raman vertices are 
of the form 



52[e,(k)+£,(k)] 



± cos (y5(k) 



dkl 



± sin (/3(k 



29^|A(k) 
dkl 



(A3) 



with the upper (lower) sign in ± corresponding to CC 
(CC). Similarly, the intra- and interband current vertices 
are 



cc,cc 



(k) 



2h 



g[£,(k)+£,(k)] 

dka 



Ji^^(k)=[J.^e(k)]* = -|-e-'^» 

. ...a^k) ,,,, 29|A(k)| 



, (A4) 



.(A5) 



The corresponding static Raman vertices 7^^''^'^(k) 
come from the effective mass theorem 



cc,cc 

f aa. 



(k) = 7^'""(k;2)± 

^ m a"%,c(k) 

~ h? dkl ■ 



2m|jg'^(k)p 



(A6) 



It is interesting also to note that for the case considered 
in Section 3 (£c(k) = -ec(k), A(k) = A^), the a = a 
vertices read as 



7Se'^^(k;2) = Tcos.^(k)7^:(k;2), 
j|XZ,cC(k) = Tcos^(k)Jr(k), 
J£^(k) = sin^(k)jr(k). 



(A7) 



where 7-(k;2) = 7-(k) = {m/h^)d^e,{k)/dkl and 
Jf (k) = {e/h) decik)/dk^ are the vertices of the A§ = 
model. 
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